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ABSTRACT 

Assuming General Relativity is correct on large-scales, Redshift-Space Distortions 
(RSDs) and the Integrated Sachs- Wolfe effect (ISW) are both sensitive to the time 
derivative of the linear growth function. We investigate the extent to which these 
probes provide complementary or redundant information when they are combined to 
constrain the evolution of the linear velocity power spectrum, often quantified by the 
function f{z)as{z), where / is the logarithmic derivative of erg with respect to (1 -I- z). 
Using a spherical Fourier-Bessel (SFB) expansion for galaxy number counts and a 
spherical harmonic expansion for the CMB anisotropy, we compute the covariance 
matrices of the signals for a large galaxy redshift survey combined with a CMB survey 
like Planck. The SFB basis allows accurate ISW estimates by avoiding the plane- 
parallel approximation, and it retains RSD information that is otherwise lost when 
projecting angular clustering onto redshift shells. It also allows straightforward calcu- 
lations of covariance with the CMB. We find that the correlation between the ISW and 
RSD signals are low since the probes are sensitive to different modes. For our default 
surveys, on large scales {k < 0.05 Mpc/h), the ISW can improve constraints on /erg 
by more than 10% compared to using RSDs alone. In the future, when precision RSD 
measurements are available on smaller scales, the cosmological constraints from ISW 
measurements will not be competitive; however, they will remain a useful consistency 
test for possible systematic contamination and alternative models of gravity. 
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1 INTRODUCTION 

Large-scale structure formation is a valuable tool for dis- 
criminating among cosmological models. When trying to 
distinguish ACDM from its alternatives, one typically char- 
acterizes large-scale structure by P(k,z), the matter power 
spectrum as a function of scale and redshift. Various late 
Universe observations, such as galajcy clustering, galaxy 
cluster counting and weak gravitational lensing, are sensi- 
tive to this function, providing snapshots of the amount of 
structure at different epochs. The dynamics of the Universe 
also produce observations that are sensitive to the rate of 
structure growth at a given epoch. Such dynamical probes 
include redshift-space distortions (RSDs) and the Integrated 
Sachs- Wolfe effect (ISW). These probes provide useful con- 
sistency checks on models of structure formation and can be 
used to distinguish models of dark energy from modifications 
to General Relativity (GR). 

RSDs are an effect seen in galaxy redshift surveys due 
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to the peculiar motions of galaxies. In the absence of pecu- 
liar motions, the redshift of each galaxy corresponds to a 
unique distance r from the observer; however, the line-of- 
sight component of a galaxy's peculiar velocity v creates an 
additional Doppler shift that distorts the observed redshift, 
thereby "moving" the galaxy to a different inferred distance 
or redstiift- space distance, s{r,v). The net effect of peculiar 
velocities is to enhance the appare nt galaxy clustering along 
the li ne-of-sight in redshift-space l|Kaiseijl 19871 : iFisher et al.l 
[1994"). On scales larger than ~10 Mpc, galaxies act as tracers 
of the bulk flow of matter and, by the continuity equation for 
a perfect fluid, the matter density at a particular position 
grows as matter flows in from surrounding regions. It fol- 
lows that peculiar velocities - and therefore RSDs - provide 
a measurement of the structure growth rate. 

At redshifts z < 0.3, the meas uremen ts from the 
2dF Galaxy Redshift Survey (2dFGRS; iCoUess et al.ii2003) 



and the Sloan Digital Sky Survey (SDSS; lYork et all 
l200d) provide the best current constr ai nts on RSD 
JPeacock et all 1200 ll: iHawkins et al.l l2003l: iPercival et al 



l2004l : iPope et al.N2004l : IZehavi et al.H2005l : lokumura et al 
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I2OO8I : ICabre fc Gaztanagal l2009al lbh. At higher redshifts, 
RSD have be en measured in the VIMOS-VLT Deep Sur- 
vey rVVDS'tjLe Fevre et al."2005|: iGarilh et al.ll200a). and 



Wiggle-Z (jPrinkwater et al.. 2QM) 

surveys (|Guzzo et alj 
l200^ : lBlake et alfcoill '). These resuhs are all fully consistent 
with the standard ACDM+GR model. RSDs also produce 
a signal when angular clustering measurements are made of 
projected d ata as R SD affect a ny redshift-d ependcnt binning 
of galaxies ijFisher ct al. 1993 : iNock et al.i,2010 ). A correc- 
tion for this has been included in previous anal yses of angu- 
lar clustering from photometric r edshift surveys l|Blake et al.l 
I2OO7I : IPadmanabhan et aLll2007l ). There is a window of op- 
portun ity where surveys such as the Dark Ene rgy Survey 
(DES; bark Energy Survev Collaboration! l2005ll, can pro- 
vide cutting edge RSD con straints Jross et al.ll201ll ') , before 
surveys such as BigB OSS (|Schlegel et al.ll2009bh and Euclid 
l|Laureiis et al.ll201ll ) extend the redshift range of spectro- 
scopic RSD measurements out to that of the photometric 
redshifts. 

The ISW occurs when CMB photons traverse evolving 
potential wells. If the gravitational potential around an over- 
density decays (grows) as a photon crosses it, the photon will 
emerge with a net blueshift (redshift), which contributes to 
the overall CMB anisotropy. In GR, potential wells in the 
late Universe are static during matter domination, reflecting 
a balance between structure growth and the background ex- 
pansion of the Universe. Once structure growth slows, due 
to e.g. the presence of dark energy, potentials decay in abso- 
lute value. The ISW is difficult to measure using the CMB 
alone but can be detected by cross-correlating the CMB 
with foreground objects th at trace the large-scale structure 
l|Crittenden fc Turok|[l99^ ). 

Initial CMB-galaxy cross-correlation measurements 
using the COBE data failed to find a signal, but this 
changed dramaticall y with the arrival of the WMAP data 
l|Bennett et al.ll2003l ). Weak correlations, at the 2 — 3a level, 
were soon seen between the CMB and numerous probes 
of large scale structure, such as th e NVSS radio galaxy 
surve y and the X - ray b ackground (|Boughn fc CrittendenI 
12004 iNolta et all 120041'). and i n optical surveys hke 



the SDSS llFosalba et al.l l2003l: IScranton et al.) 



Fosalba fc Gaztanai 



Giannantonio et al 



ga 



.2004 
20061 ). 



Padmanabhan et al 



200a 



2005 



while weaker indications 
been seen using the shallower 2MASS inf rared survey 
jAfshordi et al.ll2004 [ Passat et alllioOT: Francis fc Peacockl 



2OI0I). Combin ing these probes ( Giannantonio et al.ll20o'^ 



Ho et all 20081) raises t he signiflcance to the 4ct level (see 



also Dupe et al.l l201ll . for a comprehensive summary). 
Planck should not dramatically differ from WMAP on 
the large scales to which the ISW is most sensitive, 
but its greater resolution and frequency coverage will 
enable greater control of possible systematic contamina- 
tions from the galaxy and extra-galactic point sources. 
Improvements to ISW measurements will ultimately 
require better large scale structure data covering most 
of the sky, such as th at provided by radio surveys (e.g. 
iRaccanelli et al.l 2011) or other fr equencie s (e.g . DES 
and WISE: iDark Energv Surwv CoUabo rationf l2005l : 
IWright et alll20m 

Large redshift surveys can be cross-correlated with the 
CMB to measure the ISW effect in parallel with providing 
galaxy clustering and RSD measurements. Indeed, this com- 



bination has already been explored bv lCabre fc Gaztanagal 
(|2009ah using SDSS and WMAP data. Our goal in this pa- 
per is to understand the potential of a larger galaxy survey 
such as BOSS to constrain the growth rate of large-scale 
structure by combining the two datasets while accounting 
for the covariance between them. We will work within the 
context of GR. If one does not assume that GR holds, then 
RSD and ISW measurements contain complementary infor- 
mation: allowing for the most general metric-based models, 
ISW measurements depend on both time-like and space-like 
metric fluctuations, while RSD only depend on time-like 
fluctuations. Once Einstein's equations are assumed, both 
probes can be shown to measure the amplitude of the pecu- 
liar velocities on large scales, f{z)as{z), where 0-8(2;) is the 
normalization of the linear matter power spectrum and 



din as{z) 



(1) 



dln(l + z) 

When combining the ISW with RSDs, one naively expects 
the observables to be correlated since they both arise from 
fluctuations in gravitational potentials, but we will show 
that they are largely uncorrelated. 

For our analysis, we decompose the galaxy number 
density into discrete spherical Fourier-Bessel (SFB) modes 
(eigenfunctions of the Laplacian operator in spherical co- 
ordinates). This basis provides a natural way to relate the 
3D galaxy field to the 2D spherical harmonics of the CMB 
while avoiding approximations which may be unsuitable for 
either data set. For example, the distant observer and plane 
parallel approximations (typically applied to galaxy counts) 
poorly describe the ISW, which is only relevant for very 
large distance scales and wide angles. Projecting galaxies 
into redshift slices (typically done for ISW) discards ra- 
dial information about galaxy positions, which weakens the 
RSD signal. We could have chosen to bin very finely in 
redshift, but this would lead to a highly correlated and 
ill-conditioned covariance matrices. Previous works have 
also found the SFB basis useful fo r analysing; various 
probes of large-scale-structure (Fish er et al. 19951: 
Schmoldt et al. 1999; Heavens 2003 ; ICastro et al.l 
2005; Erdoedu et al. 2006a|Jb|; lAbramo et al.l I2OI0I : 
Rassat &: Refregier,,2011i ). 



Outline 

We begin in [J2]with a brief summary of the survey assump- 
tions and the cosmological model on which we base our sub- 
sequent calculations. In §31 we review the SFB expansion 
and apply it to the ISW signal and to the galaxy number 
density in redshift-space. In 33] we compute the covariances 
and correlations of these signals. In §S1 we forecast the abil- 
ity of RSDs and the ISW to constrain the growth rate when 
combined. We conclude in ^ The appendices discuss issues 
related to the sampling of discrete fc-modes in the SFB for- 
malism. 



2 DESCRIPTION OF FIDUCIAL SURVEY 

We consider a spectroscopic redshift survey which measures 
the redshift z and sky position Q for a large sample of galax- 
ies. For simplicity, the survey is assumed to cover the full 
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solid angle on the sky and have a mean observed comov- 
ing galaxy density p^^^{z) that is constant over the redshift 
range, (0, Zmax). We assume that galaxy redshifts and densi- 
ties are similar to those expected for the BOSS survey, and 
take 2max = 0.6 and p^''\z) = p*^"' = 3 x 10""* / Mpc^ 
with Ho — h X 100 km/s/Mpc. Our fiducial survey sub- 
sequently has 4 times the volume of BOSS, which will only 
survey 1/4 of the sky. The galaxy sample is assumed to obey 
a scale-independent linear bias with respect to dark matter 
clustering: 



(2) 



where we ignore the gauge-dependent effects of the bias def- 
inition in GR (these introduce cor r ections for horizon-size d 
scales; see IChallinor fc Lewisll201ll : iBonvin fc Durredl201lD . 
We also assume constant galaxy clustering (CGC), i.e. that 
the observed power spectrum of galaxy density P^^'(fc) is 
constant in z. Physically, CGC means that the galaxy bias 
is related to the linear growth rate by 



b{z) 



m 

Giz) 



(3) 



where matter density modes grow by G{z) = S{k; z)/S{k; z = 
0) in linear theory. Linear bias and CGC a re suitable ap- 
proxim ations for BOSS, for which b{0) = 1.7 (jSchlegel et al.l 
l2009al ). We restrict our analysis to galaxy density modes 
with k ^ 0.1/i/Mpc; for our purposes, these modes are de- 
scribed sufficiently accurately by linear theory. 

For Planck, we assume all-sky, cosmic-variance limited 
measurements of the CMB temperature anisotropy. We re- 
strict our analysis to large angular scales (>1 deg) for which 
details of the beam will be irrelevant. 

Our fiducial cosmological model is a flat ACDM model 
with the following parameters: 



Q,m f^b h n as 
0.27 0.045 0.72 1 0.75 

The parameters ilm and Qt are the total matter and bary- 
onic matter densities as fractions of the critical density; n 
is the tilt of the primordial spectrum of scalar fluctuations; 
and erg is the normalization of the matter power spectrum 
today, expressed as the RMS of linear density fluctuations 
smoothed on scales of 8 Mpc//i. All figures and calculations 
use these parameters unless otherwise noted. We work in 
the Newtonian gauge, for which the perturbed Friedman- 
Robertson- Walker metric is described by 



ds 



a[ -(1 + 2ip)d7j'^ + (1 + 2(l>){dr + rdVl) 



(4) 



where r) is the conformal time, (j) and i/; are the grav- 
itational potentials, and a = (1-1- z)"'^ is the cosmic 
scale factor. We calculate the linear matter power spec- 
tr um (without bar y on w iggles) using the fitting formulae 
of lEisenstein fc Hul (|l999l ). 



2D CMB anisotropy. This approach allows a straightforward 
calculation of covariance matrices for the two observables, 
as we show in fHl 



3.1 Review of the SFB expansion 

The SFB expansion allows us to express any scalar func- 
tion in terms of the eigenfunctions of the Laplacian opera- 
tor in spherical coordinates. Galaxy number density can be 
expanded thusly: 



(5) 



!=0 m — — l 



where r = (r, 0) is the comoving distance vector, ji is a 
spherical Bessel function, Yim is a spherical harmonic, and 
cin is a normalization constant defined below. For a given Z, 
the ki n and Cin will depend on a choice of boundary condi- 
tions. IHeavens fc Tavloil (Il995l ) proposed using the discrete 
wavenumbers that extremize the spherical Bessel functions 
ji{kinr) on the survey boundary (Neumann boundary con- 
ditions) : 



dr 



ii(kinr) 



0, 



(6) 



where rmax = '"(zmax). This choice avoids redshift-space dis- 
tortions on the survey boundarjQ, i.e. the survey boundary is 
a sphere of radius rmax in both real-space and redshift-space. 
With this choice of boundary conditions, Sturm-Liouville 
theory guarantees that the j;(fc;„r) form a complete basis 
for the r adial part of the expansion on the finite interval 
(0,rmax) (|Wang et al.l l200g). 

For our full-sky survey, the range of I and n required is 
determined by the radial distribution of galaxies and the de- 
sired 3-D resolution. The / value determines the wavenum- 
ber perpendicular to the line of sight, while the fc;„ indi- 
cate the wavenumber along the line of sight. For a fixed 
3-D wavenumber, fcmax, higher I values have fewer line of 
sight modes. For each l^2,we order the discrete wavenum- 
bers such that kin monotonically increases with n, and we 
define the smallest kin to have n — 1. For / = 2, we 
find that we need all n up to n ~ rmaxfcmax/7r to ac- 
count for all k2n < fcmax. Larger I require fewer n, and 
the largest included / may have only one mode {n — 1) 
below the maximum k. In our case, in order to include all 
kii < femax, we find that we need / up to / « rmaxfcmax. The 
smallest k (longest wavelength) included in the analysis is 
fc2i — 3.342/rmax; for our fiducial survey with Zmax = 0.6 
and r-niax = 1563 Mpc//i, we have k2i = 0.00214 fe/Mpc. 

The spherical Bessel functions obey the following or- 
thogonality relation: 



dr ji{kinr)ji{kin'r)r^ = c,„^5w 



where 5 is the Kronecker delta, and 



(7) 



3 THE SPHERICAL FOURIER-BESSEL 

EXPANSION APPLIED TO LARGE-SCALE 
STRUCTURE 

In this section, we use a SFB expansion to relate the matter 
power spectrum to both the 3D galaxy density field and the 



^ The boundary conditions restrict our ability to reconstruct cer- 
tain real-space features of the density field; however, we are only 
interested in the power spectrum of modes that are sampled by 
the survey, not accurate reconstructions of the density and veloc- 
ity fields. See Appendix lAl for further discussion. 
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jl (^in^max) 



The Spherical harmonics are also orthogonal: 



Therefore, the SFB coefficients are given by 



3.2 



(8) 



(9) 



(10) 



Galaxy Density Fluctuations in 
Redshift-Space 



In a flat Friedmann-Robertson-Walker (FRW) metric, red- 
shift and comoving distance are related by ^ = H{z), where 
H{z) = a /a IS the Hubble parameter. This relation is applied 
to measured redshifts to obtain a redshift-distance, s, which 



includes distortions from peculiar velocities. Percival et al 



2004 , hereafter P04) applied the work of lHeavens fc Tavlor 



19951 , hereafter HT95) to the 2dFGRS, providing details 



of a practical implementation of the SFB analysis to galaxy 
counts. Here, we reproduce the calculation described by P04, 
simplified for an all-sky survey with CGC, no weighting, and 
no non-linear RSDs on small scales (Fingers-of-God). In the 
notation of P04, Greek subscripts refer to a triplet of indices 
v = (l^m^n^) and>(r) = ji^{ki^„^r). 

The galaxy overdensity in redshift space, D{s) = 
p^^^{s) — p^^'(s), can be transformed into its SFB coeffi- 
cients, 

Dl=c, j f "d^s D{3)j.{s)Y:{e) , (11) 

J iir J 

and expressed as linear combinations of the real-space mat- 
ter overdensity coefficients, 5^, at a = 0: 



i?; = ^[6(o) 



,/(0)Km]5^ 



(12) 



There is a usual clustering term 

X ^""dr r'Mr)j,ir)r''\r)G{z)^^^ 
and a RSD term 



(13) 



0;. .;..£)„ 



(14) 



where z = z{r) is evaluated on our past light-cone and 
j' = dj /dr. Note that for an all-sky survey with no masking, 

is symmetric but VL;j is not: f{Q)Vi,p, is the contribu- 
tion to the observed galaxy density mode Dl from velocities 
source d by matter m ode 5^^. We will use the approxima- 
tion of iLindeil l|2005h . who showed that f{z) r; Q,„,{zf-'^^ 
for ACDM. (other approximations have been derived by e 
|Peebles|[l980l : iLahav et al.lll99ll : IWang fc SteinhardtllT99; 

Since our simplified $ and V matrices are independent 
of m and block-diagonal in I and m, the are more simply 



expressed as 



cinCi„t I dr r ji{rkin)ji{rki„, 



V. 



X ■p^''\r)G(z) 



k? , 

Ln' 



m 

dr r"^ Ji{rkin)j[{rkir. 



xr^\r)G{z)ll^. 



(15) 



(16) 



(17) 



Using Eq. (O, the CGC condition of Eq. ((3|, and the orthog- 
onality of the spherical Bessel functions, we find that "l>5j,i/ 
is simply proportional to the identity matrix when p^^'(2) 
is constant: 



dr r jt{rki„)ji{rki„>) 



— eal r —2 cK i — e-al cK 



(18) 



Sections of the RSD matrix V^„, are illustrated by 
solid lines in Figure [T] (assuming CGC). Even in our ide- 
alized survey, the matrices exhibit significant mode mixing 
at low I. Mathematically, the negative off-diagonal elements 
of V^„, arise because, unlike the spherical Bessel functions 
in Eq. pSI) . the Bessel function derivatives in Eq. (|17|) do 
not form a perfectly orthogonal basis. This result contrasts 
with the familiar Kaiser result (|Kaisei|[l987l ). which has no 
mode mixing: 



D^k) =S{k) [6+(fc-r)V] 



(19) 



Kaiser approximates the radial distortions as being along a 
particular Cartesian axis, which allows a plane-wave expan- 
sion. Whereas the derivative of a single plane-wave is an- 
other plane-wave with the same frequency, spherical Bessel 
functions do not enjoy this property. Therefore radial mode- 
mixing is an inevitable geometric effect of RSDs on large 
angu lar scales (Zaroubi fc Hoffman 199^: lHeave"ns^ fc Tavloil 
Il995h . For the correlation function and its Legendre mul- 
tipole expansions, the full wide angle effects including this 
mode mixing can be calculated using a tri-polar s p herica l 
harmoni c s exp an sion as described by ISzalav et al.l (|l998f ): 
Szapudil ||2004|): IPapai fc Szapudil l|2008ll . and tested by 



(2010|). Although we restrict our analy- 



sis to < 0.1/i/Mpc, the observable redshift-space modes, 
Dimm with kin just below this cutoff will have contribu- 
tions from real-space modes, Simn' with fc;„/ just above 
the cutoff. Therefore, computations must extend beyond 
the cutoff scale to ensure convergence. We find that includ- 
ing modes up to = 0.15ft/Mpc is more than sufficient 
for our purposes and that we can safely ignore mixing be- 
tween modes with Afc > 0.05/i/Mpc; however, surveys with 
a non-trivial p^^'(z) or selection function will exhibit greater 
mode-mixing. HT95 describe how to attenuate mode- mixing 
by optimally weighting the data. 
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Figure 1. Sections of the V^^, and U'^^, matrices for 1 = 2 (top) 
and Z = 20 (bottom); respectively, these give the contributions to 
the RSD and ISW signals from a real-space matter density mode, 
^Imn' ■ For clarity, V^^, for n=20 has been offset by -0.3, and all 
matrices have been scaled to fit on a common vertical axis (see 
legend). The central spikes of the V^^/ occur where n' = n. Note 
that for a given I, the lowest available wavenumber is fcji- 



3.3 The Integrated Sachs- Wolfe Effect 

The ISW c ontribution to the 2D CMB temperature 
anisotropy is (|Sachs fc Wolfdll967l ') 




T 



Arj e 

' va 

rvi,s 
2 d-qe 



d 



drj 

drj 



(20) 



where rj is conformal time, t^o is the present time, t^ls is the 
time at the last scattering surface. The last equality follows 
from assuming GR and the absence of anisotropic stress, so 
that tp = —(f). The integration is done along a line-of-sight 
trajectory r = (r, G) where r = r{ri) is the comoving distance 
to an observed CMB photon at time r;. Although the integral 
is along our past lightcone, note that the time derivative ^ 
is simply taken with respect to the global time coordinate, i]. 
The CMB scattering function r(r;) is much less than unity, 
so we will neglect the factor hereafter. Thus, we take 



AT'' 



T 



-{9) 



vo 



drj 



'^'^ 5 

or 



(21) 



where we have changed variables to match the previous sec- 
tion. The z argument above reminds us that (j) is being eval- 
uated along our past light-cone at z = z{r). 

Although the ISW arises from all potentials between 
the observer and the last-scattering surface, surveys are pri- 
marily sensitive to potentials within the survey boundary, 
which is a sphere of radius r-max = r{rjma.^). Therefore we 
write 



-(0) 



d^MZIil + sTdrMlIi), (22) 



dr 



dr 



and ignore the second term on the right since it has only 
a relatively small correlation with the galaxies in the sur- 
vey. We label the remaining term with a "CC" for "cross- 
correlate" : 



AT^ 



T 



-{0) 



dr 



(23) 



Expanding the potential into SFB modes, we find 



AT' 



iCC 



T 



-{0) 



10 



or 



E 

Imn 



{z)ji{kinr)Yim{0) 
d 



= 2^ci^Yi^{0) j dr ji{ki„r)—(f>imn(z){2'i) 

Imn 







where the ki„ and ci„ are (conveniently) the same as those 
for the galaxy survey. The z argument remains in (j>imn{z) 
to account for linear growth. 

Einstein's equations allow us to convert potential to 
density using Poisson's equation0 

\7^(f){r) = -47vG5p{r) , (25) 
which simplifies substantially in the SFB basis: 

(l>lmn{z) = -2-47rG Spimniz) . (26) 

The scale factor a = 1/(1 + z) accounts for the metric in an 
expanding Universe (see e.g. lDodelson„2003i ) . We can rewrite 
this as 



Plmn\Z) — — lim-nQ- 



Ln 



(27) 



where 5 is the present comoving matter overdensity and G{z) 
is the growth factor (not to be confused with the gravita- 
tional constant, which been absorbed into H^). Plugging 
Eq. dm into Eq. JgU and defining 



CC 



A T-iCC 

d'0 -^{0)YUO) 



we find that it is a linear combination of the Simn: 

G{z) 



CC 



SrimJ/o ^5;m.n / drji{kinr)-^ 

' ^In Jo 



(28) 



.(29) 



To reiterate, af^ is not the full ISW contribution to the 
CMB anisotropy - it's just the contribution from density 
modes in our analysis, which are restricted by the galaxy 
survey geometry and our boundary conditions. 
We now rewrite af^ in terms of f{z): 

CC _ on \ ^ r 



dr ji{kir,r)H{z)—[{l + z)G{z)] 



^ Strictly speaking this equation depends on the gauge where the 
density fluctuations are defined, but on sufficiently small scales 
all gauges are equivalent. For the bulk of the ISW signal-to-noise, 
this is a good approximatio n, but the low I results can be affected 
llYodl2009l: IYoo et aljlioogh . 
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Figure 2. Functions related to tiie linear growtii of large-scale 
structure. RSDs are sensitive to f(z)G{z) while the ISW is sen- 
sitive to [1 - f{z)]Giz). 



X / Arji{ku,r)H{z)G{z)[l-f{z)] (30) 



where we have used / — 
convenience, we define 



dln(l + z) dr 



'^In, Jo 

SO that 



dr ji {ki„r)H{z)G{z) [1 - f{z)] (31) 



(32) 



Sections of the t/^ matrix are shown in Figure [T] 

3.4 Comparison of Sensitivity to Large-Scale 
Structure Growth 

Notice by Eq. ^ and Eq. ^ that RSDs and the ISW are 
sensitive to complementary functions of the linear growth 
and growth rate: fG{z) and (1 — f)G(z), respectively. Both 
functions can be expressed in terms of fas{z), since 



asiz) = cr8(0)G(z) 



(33) 



There is a relative H{z) factor between the RSD and ISW 
dependences; however, these data will not be a primary 
probe of that function. As shown by the dashed curves in 
Figure O the ISW signal is strongest at late times when the 
Universe is becoming A dominated and gravitational poten- 
tials are decaying most rapidly. The RSD signal is strongest 
at intermediate redshifts but maintains a relatively constant 
signal over the range plotted. This reflects a balance be- 
tween f{z), which drops as the Universe leaves matter dom- 
ination, and G{z), which continues to increase. The real- 
space galaxy clustering, described by Eq. ([T6|l . is sensitive 
to bag{z), which is constant in the case of CGC. 

Since both the RSD and ISW signals can be expressed 
as sums over the matter density modes, Simn, we generally 
expect some correlation between them. As we show in the 
next section, the degree of correlation is determined by a 
contraction of the and V^„, matrices. Due to their rela- 
tively orthogonal shapes (shown in Figure[T]), the correlation 



turns out to be low, which makes the ISW and RSDs inde- 
pendent probes of the structure growth rate despite being 
based on closely related underlying physics. 



4 COVARIANCES AND CORRELATIONS 

Having written down expressions for the redshift-space 
galaxy density modes Df^^ and the ISW contribution to 
the CMB anisotropy aj^, we will now calculate their co- 
variances - measurable quantities that we can predict. Since 
Di^„ and af^ are both linear combinations of the real- 
space matter density modes, Simn, the covariance between 
any two observable modes will be a linear combination of 
(5;mn<5;*m/,i/), where angle-brackets denote an ensemble av- 
erage. Fortunately, in linear theory, 



(34) 



Hence, all covariances will simplify substantially and be ex- 
pressible in terms of P(k). By the symmetry of the <I>, V, 
and U matrices, the covariances between observable modes 
will vanish unless I = /' and m = m' , so we can restrict our 
attention to these cases. 

Using Eq. ((TSJ and Eq. p4p , the covariances of the den- 
sity modes are 



Imn Imn / 



j2p{kt„„)[b{o) 



+ /(o) vi, 



X [6(0) $l„„+f{0)V:,,„„] 



(35) 



Note that the order of the subscripts of V^^, is impor- 
tant since that matrix is asymmetric. The measured modes, 
Dimn: will include shot noise; for a constant-density/all-sky 
survey, the total covariance will be 

{Di^„Di^j^i'^ = {Di^„Di^„i) + A,j„/ . (36) 

In our simplified case of an all-sky survey with no weighting 
or masking, the shot noise term is given by A'^^, = . 

We cannot measure the ISW signal, af^, in isolation: 
we must correlate the galaxy density with the total mea- 
sured CMB anisotropy, which includes contributions from 
the early Universe and detector noise. That is, we can only 
measure (^aimDtmn)- However, the predominant galaxy- 
CMB correlations are sourced by gravitational potentials 
within the survey boundary, so we can safely ignore CMB 
anisotropies originating at higher z and on the surface of 
last-scattering (chance correlations with these sources will 
be accounted for in our covariance matrix). Furthermore, we 
do not expect correlations between the noise of the two types 
of observables. Hence, we let l^aimDtmn) = {'^?mDi*nn}- 
Combining Eq. ((SH , Eq. JTHJ , and Eq. JMl , we find that 



(ap™ A™) = J2 PMHo) $L' + /(o) vL']uL> . (37) 

n' 

Note that the VU term is the "efltect of RSDs on the ISW 
signal," which we discuss below. The covariances of the af^^ 
are 



(38) 



Again, these aren't directly observable, but they contribute 
to the total CMB covariance, (^\aim\'^')- 
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Figure 3. Top: Pearson coefficients for CMB-galaxy cross- 
correlations. The solid line includes the effects of RSDs while the 
dotted line ignores them. The are calculated using the ISW 
signal alone, and the Df neglect shot noise, therefore, these 
Pearson coefficients are theoretical upper limits for the signals 
alone. Bottom: Pearson coefficient computed using the RSD term 
only (6 = 0). The curve oscillations are not due to numerical noise 
- they arise from the , matrix. 
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Figure 4. The change in CMB-galaxy cross-correlations (top) 
and galaxy auro-correlations (bottom) due to RSDs. is the 

same as D"^^ defined in Eq. II15I I but with the RSD term re- 
moved. The change in the cross-correlation is typically small, ex- 
cept on the largest angular scales where it can be significantly 
negative. The curve oscillations are not due to numerical noise - 
they arise from the , matrix. 



It is interesting to look at how the presence of RSDs 
afTects the correlation between the galaxy density and the 
CMB anisotropy. What we need are the Pearson correlation 
coefficients, defined as 



RxY 



(XY) 



(39) 



for two random variables, X and Y . Two modes are fully cor- 
related (contain completely redundant information) when 
their Pearson coefficient equals ±1, and they are completely 
independent if the Pearson coefficient is 0. We plot CMB- 
galaxy Pearson coefficients in the upper panel of Figure |3l 
Note that the and af^ plotted do not include noise or 

contributions from the non-ISW CMB. These Pearson coef- 
ficients are therefore a theoretical upper limit - the measur- 
able quantities and aim can only be less correlated. 

Figure |3] shows that the inclusion of RSDs reduces the 
CMB-galaxy Pearson coefficients by a small amount for our 
modes of interest. We highlight two reasons for this decrease 
in Figure [l] On small angular scales, it is primarily because 
RSDs increase the observed galaxy auto-correlations, while 
the CMB cross-correlations are left relatively unchanged by 
the RSD term in . Therefore, the effect on the Pearson 
coefficient calculation is that RSDs increase the denomina- 
tor more than the numerator. On the largest angular scales, 
CMB-galaxy correlations are also reduced because the ISW 
can correlate negatively with the RSD part of galaxy cluster- 
ing. To see why this happens, note the negative side-lobes of 
the V^^i matrix in Figure [T] which are due to mode- mixing. 
The Un matrix for the ISW is always positive, and when con- 
tracted with V^^i , the negative contributions to the sum can 
dominate. Hence the reduction in CMB-galaxy cross corre- 
lations is partially a geometric effect of the RSDs on large 
angular scales. 

To further demonstrate the lack of correlation between 
the ISW and RSD signals, we can compute the Pearson co- 
efficients for CMB-galaxy correlations using only the "dis- 
tortions" - the observed galaxy density fluctuations arising 
only from the RSD term in Eq. (|12|l . This computation is 
done by setting fe(0) = 0. Thus we are imagining a galaxy 
density field with no real-space clustering but with observed 
redshift-space clustering due to peculiar velocities arising 
from the gravitational potentials. It is clearly an unphysi- 
cal situation; nevertheless, this Pearson coefficient allows us 
to understand where information about structure growth is 
coming from. We want to discard the real-space galaxy clus- 
tering term, which has a very high signal-to-noise, but which 
is insensitive to /erg (2). We see in the lower panel of Figure 
[3] that correlations between the ISW and RSDs do exist via 
the "VC" term in Eq. ((37]), but they are at most 20% on 
the largest angular and radial scales. As before, this is a the- 
oretical upper limit excluding noise and the early Universe 
CMB anisotropy, which can only reduce the correlations. 
This low correlation, even in the absence of noise, occurs 
because the ISW and RSD signals depend on rather orthog- 
onal combinations of the matter density modes under GR. 
Since RSDs directly probe the time-like potential tp while 
the ISW is sourced by t/j — 0, one might have hoped to use 
the two signals to test GR by measuring anisotropic stress 
on a mode-by-mode basis. The low correlation between the 
signals implies that such an approach would not be feasible, 
but it does not prevent statistical tests of GR. 
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We further find that the presence of RSDs does not af- 
fect the signal-to-noise ratio for an ISW detection. To com- 
pute the signal-to-noise, we first define Cf^ = l^Dimn'hm)- 
The covariance of Cf^ is 

Cov(C"'^ (joT-j ^ ^In '^L' + (kimP) {PlmnDlmn') ^^^^ 

and the total signal-to-noise is given by 
(I) ' = E E CfJCov-((7f:, CfJ) (41) 

I nn' 

For our fiducial galaxy survey and Planck, we find that the 
total ISW signal-to-noise ratio is 4.25, and this val ue is neg- 
ligibly affected by including or excluding RSDs. lAfshorda 
('2004) calculates that an ideal survey extending to a red- 
shift of 2 to 3 could detect the ISW at about 7.5cr. The 
limiting factor for our signal-to-noise ratio is the cutoff of 
the galaxy survey at Zmax = 0.6, and our result is consistent 
with Afshordi's findings at the lower redshift. 




0.01 1 . . , , . . , . . ■ ■ 'I 
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5 FORECAST FOR COMBINED 

MEASUREMENTS OF THE GROWTH RATE 

In this section, we forecast our ability to use galaxy and 
CMB maps to constrain /ug - the normalization of the linear 
velocity power spectrum - as a function of redshift. Under 
GR, f{z)iJi{z) « ijs,{z)^m{z)''''^\ where 7(2:) is nearly con- 
stant and weakly dependent on the dark energy equation of 
state (lLinderll2005l:|Polarski fc Gannouiill2008lJPeebleslll980l : 
iLahav et al.lll99ll : IWang fc Steinhardtlll998l '). Instead of fo- 
cusing on dark matter and dark energy parameters, we fore- 
cast direct constraints on /erg itself in a few redshift bins. 
Recall from H3. 41 that /ug is mainly probed by the ISW and 
by the RSD term in the redshift-space galaxy density. The 
real-space galaxy clustering also constrains h{z)(7s[z)\ how- 
ever we are less interested in this function, which is a con- 
stant in the case of CGC. 

In general, for a set of cosmological parameters that 
we wish to constrain, the Fisher matrix of the parameters is 



Fc,p = -Tr 



(42) 



where C is the full covariance matrix for all observables. The 
forecasted errors on the pa are (Tfix(pa) = Faa^'^ if all other 
parameters are held fixed or (7marg(Pci) = {F~^)l/a if we 
marginalize over the other parameters. 

Our observables are the aim, the spherical harmonic 
coefficients of the CMB anisotropy, and the -Df„„ , the SFB 
coefficients of the galaxy number density. We combine these 
into a joint data set, Aimn'- 



air. 
DL 



(n = 0) 
(n > 0) 



(43) 



The joint covariance matrices are Cuirnm'nn' = 
<^AimnA;*^,jj,^. Due to the symmetry of our fiducial 
survey and the linearity of our cosmological model, the 
C will be independent of m and block-diagonal in m and 
I. Therefore, we need only calculate the much simpler 



Figure 5. Top panel: Forecasted error bars on pa, the offsets to 
/(Tg, as a function of maximum k. The solid, dotted and dashed 
curves are for the low, medium and high redshift bins, respec- 
tively. For each pair of curves, the upper curve uses galaxies alone 
while the lower curve includes CMB-galaxy cross-correlations. 
The constraints on each z bin have been marginalized over the 
other two bins and the normalization 6(Tg(2 = 0). Bottom panel: 
For each pair of curves in the top panel, the ratio (upper over 
lower) minus one. 



submatrices, C'j„, = (Aim„A*^,j,), and we find that 



I 

= 5E(2^ + i) 



E(^ 



1 ^^n'n" /r>l\-l ^^n"'r 



nn' 



' dp^ 'r^"n'" 



(44) 



The Pa for which we forecast constraints will be offsets 
to fas - relative to ACDM - in three redshift bins. We 
divide our redshift range into three bins of equal width in 
log(l -I- z). The bin boundaries are zo = 0, zi = 0.170, 22 = 
0.368,23 — 0.6. We start with a fiducial /ag(z) computed 
using ACDM and allow this function to have a piece-wise 
constant offset in each bin: 

/ag(z) = /erg (z; ACDM) {zc-i < z ^ z^) (45) 

for a=l,2, or 3. We do not specify fas{z) beyond the survey 
boundary {z > Zmax), but it is constrained by the normaliza- 
tion of the CMB. When we modify /c7g(z), we then compute 
the linear growth G{z) by integrating fG = gf^, normal- 
izing so that (Tg(0) remains fixed. Additionally, we include 
the normalization of the present galaxy power spectrum as 
a nuisance parameter: po = 6crg(0). Although RSDs and the 
ISW have some sensitivity to other cosmological parame- 
ters, such parameters are left fixed in our forecast since we 
expect them to be tightly constrained by other observables. 
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For instance, the ISW is sensitive to Q,mh? , but this param- 
eter is well-measured by the first acoustic peak of the CMB. 
We compute using two-sided derivatives with steps of 
Apa =0.01 or 0.005 and find little difference in the results. 

Figure [5] shows our forecasted constraints on the pa, 
with and without CMB information, as a function of our 
cutoff scale for the observed galaxy density, fcmax- The er- 
rors have been marginalized over the other parameters, in- 
cluding the normalization h(jg,{Q). The bottom panel of Fig- 
ure [5] shows the extent to which using the ISW has im- 
proved the constraints on fas{z) relative to using RSDs 
alone. The plot shows that for k < 0.05/i/Mpc, ISW mea- 
surements would improve constraints in each redshift bin by 
more than 10%. When RSD measurements are available on 
smaller scales, the ISW does not significantly improve con- 
straints on a scale-independent /erg (2). Our constraints for 
kTa a.yi = O.lfe/Mpc ar e in good agreement with the results 
of (| White et al.|[2009l '). who use the flat-sky and distant ob- 
server approximations. To test the sensitivity of our results 
to boundary conditions, we tried recomputing the fas con- 
straints using different values of Zmax (0.6, 0.8 and 1.0) while 
keeping the definition for pa in Eq. (|45|) . We find that our 
results in Figure \S\ are reasonably robust, and summarize 
them in SjX] 



6 CONCLUSIONS 

The integrated Sachs- Wolfe effect and redshift-space distor- 
tions are unique cosmological probes in that they are sen- 
sitive to the evolution of large-scale structure rather than 
simply the amount of structure at a given epoch. In gen- 
eral, they probe different combinations of the linear gravi- 
tational potentials, (p and tp; however, under GR, they both 
can be shown to measure fas{z), the amplitude of the pecu- 
liar velocity power spectrum. We investigated whether the 
two probes were complementary or redundant when com- 
bined to constrain fas{z) using data from the Planck CMB 
experiment and a large galaxy survey with a redshift distri- 
bution similar to that expected for BOSS. Our analysis used 
the spherical Fourier-Bessel expansion in order to account 
for important large-angle effects and to avoid discarding in- 
formation via redshift-binning. The SFB basis also provides 
an original insight into the effects of RSDs on CMB-galaxy 
cross-correlations, and it allows us to make a mode-by-mode 
comparison of the correlation between RSD and ISW mea- 
surements. 

We find that the ISW and RSDs are mostly independent 
(uncorrelated) observables, being sensitive to mostly orthog- 
onal combinations of the matter density field. In general, 
RSD depend on radial modes, while the ISW is more sen- 
sitive to angular fiuctuations in the galaxy density. We also 
find that the ISW, measured through CMB-galaxy cross- 
correlations, improves fas{z) constraints from RSDs by only 
10% even when the analysis is restricted to large physical 
scales {k > 0.05/i/Mpc). Thus, when future precision mea- 
surements of RSDs are available, the ISW will be more valu- 
able as a probe of non-standard GR models and as a test 
of survey systematics and less valuable as a way to measure 
cosmological parameters in GR. 




0.001 0.010 0.100 

fc,„ (hMpc) 

Figure Al. Histogram showing how P{k) is sampled by the 
discrete ki„ of the SFB expansion. The dark and light solid 
lines show the number of ki„ in bins of Afc = O.OOlh/Mpc for 
^max = 0.6 and Zmax = 1, respectively. Both assume Neumann 
boundary conditions. The linear matter power spectrum is shown 
by the dashed line. 
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APPENDIX A: i^-SPACE SAMPLING 

The discrete fc;„ needed for our SFB analysis are determined 
by the boundary conditions. Figure lAll illustrates the sam- 
pling in fc-space corresponding to the Neumann boundary 
conditions in Eq. Changing ^max changes the values of 
the kin, the number of kin less than fcmax, and subsequently 
how P{k) is sampled. In particular, the minimum wavenum- 
ber is given by ^21 = 3.34/rniax, which in our case is 0.00214 
/i/Mpc for Zniax=0.6 and 0.00142 /i/Mpc for z^ax^l. Since 
our Fisher matrix analysis is sensitive to P{k), it will be 
somewhat sensitive to our boundary conditions. 

We tested the sensitivity of the fas constraints in Fig- 
ure [S] to boundary conditions by varying the Zmax which 
defines the SFB expansion. We choose to compare the un- 
marginalized errors for different boundary conditions since 
increasing the survey size significantly improves constraints 
on bas, which is partially degenerate with the fas con- 
straints. For Zmax = 0.8 and fcmax = 0.1, we find that the 
unmarginalized error bars on fas change negligibly in the 
two lower- 2; bins, and they degrade by less than 2% for the 
highest z bin. For Zmax = 0.8 and femax ~ 0.02, the errors 
degrade by at most 10% for the high-z bin (with ISW) . The 
larger effect for the lower fcmax reflects the fact that the lower 
k are more sparsely sampled; therefore changes to the sam- 
pling of -P(fe) have a greater impact on constraints relying 
on the lower k. Pushing the boundary further to Zmax ~ 1 
results in negligible changes from Zmax = 0.8. 
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space. The spherical harmonic transform of D2 (0) is 
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Figure Bl. The effects of boundary conditions on computations 
of CMB-galaxy angular cross-correlations (top) and galaxy auto- 
correlations (bottom) in redshift slices. The denominators, Cf"^ 
and Cj^^, are computed by summing coefficient correlations in 
Eq. IIB5II and Eq. IIB6II . with the boundary set at ^max = 1-0. 
AC; is Ci{zma.x = 0.8) — C;(2max = 1.0). The redshift slices are 
Gaussians centered at zq with cr^ = 0.02. 



We could have chosen to make Zmax '3> 1 so as to 
have the finest possible sampling in fc-space; however, doing 
so adds considerable computational difficulty without much 
benefit. In addition to increasing the number of wavenum- 
bers, our galaxy survey would require a selection function 
that drops the galaxy density to zero above z ~ 1, resulting 
in very significant mode-mixing in the $ and V matrices in 
Eq. dSI- 



APPENDIX B: ANGULAR CORRELATIONS IN 
REDSHIFT SLICES 

A more familiar method of computing correlations of galaxy 
number counts is to bin galaxies into redshift slices. The 
2D projected galaxy density in each slice can then be auto- 
correlated or cross-correlated with other slices or the CMB 
temperature. In this section, we show how to convert our 3D 
SFB results into the more common 2D representation. 

Start with the redshift-space galaxy density contrast 
expanded into SFB modes: 



D{s) = J2 '^'"■O'" 



with 



<f9YiUe)D2{e) 



As w{s)ji{kins) 



(B3) 



(B4) 



Correlating these projected density modes with the CMB 
temperature shows that the 2D/3D CMB-galaxy correla- 
tions can be summed accordingly to obtain their fully 2D 
angular correlation in the slice: 



= (^04rnDim) — ^ ^ Win (llrn-Dfrnn) 



(B5) 



As before, the symmetry of our fiducial survey ensures inde- 
pendence from m and delta functions in I and m. Similarly, 
the autocorrelation of galaxy clustering can be written 

Cf« = (DLCf™> = Yl W^inW/,;, (i3Ln A" n' > • (B6) 



Computing Cf ® and Cf using the above sums is a bit 
challenging since since their accuracy depends on the bound- 
ary conditions (see SjAj and on the maximum k. The expres- 
sions should converge for large enough Smax and k, but the 
computation becomes intensive due to mode-mixing and the 
increasing number of discrete modes. More manageable ex- 
pressions for computing galaxy correla tions in redshift slices 
(with RSDs) are given by Padmanab han et al.l (|2007l ) and 
iRassatI (|2009l ). 

To demonstrate the limitations of our approach, we 
compute Cf ^ and Cf for several redshift slices and bound- 
ary conditions. We use Gaussian slices with (z) = zo and 
(Jz = 0.02 so that the window function w oc exp[— i(2; — 
zo)^/c"z]- Figure iBll shows how our angular correlation cal- 
culations change when we change boundary conditions from 
Zmax = 0.8 to Zmax = TO. We sec that both the galaxy count 
auto-correlations and the CMB-galaxy cross-correlations are 
affected by several percent on the largest angular scales. This 
is because the characteristic scale of a slice, k ~ l/r{zo), lies 
in a poorly sampled region of the ki„ for small I (see Fig- 
ure [XT|. For example, with I = 2 and zq = 0.6, we have 
l/r{zo) = 0.00128 /i/Mpc. We reiterate that these boundary 
effects have only a small impact on our main result, which 
is the forecast in 



,rjl{kl„s)Ylm(0) 



(Bl) 



The projected density contrast in a redshift slice is 

D2(e) = /" ds w{s)D{s) 
Jo 



(B2) 



where ■w{s) is a radial window function defined in redshift- 
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